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The interlayer exchange couphng (lEC) of two local moment ferromagnetic layers separated 
by a non- magnetic spacer layer (M/N/M multilayer) is studied using the modified RKKY method 
along with the s-f model. The TEC exhibits oscillatory behaviour with respect to the spacer layer 
thickness and it oscillates between ferro- and antiferromagnetic configurations. The conventional 
RKKY method is also used to obtain the lEC and the results are compared with those obtained from 
the modified RKKY method which incorporates the electron correlation effects. We find significant 
correlation effects on the lEC and in fact the correlations alter the nature and magnitude of the 
magnetic coupling. Hence this study points out the importance of the inclusion of correlation effects 
in the understanding of the lEC in multilayer systems with local moment ferromagnetic sub-layers 
and to offer satisfactory explanation for the experimental results. 



I. INTRODUCTION 

The indirect exchange interaction between two fer- 
romagnetic layers separated by a non-magnetic spacer 
layer exhibits oscillatory behaviour with respect to the 
spacer layer thickness and this has been observed ex- 
perimentally in many magnetic multilayer systems [1-6] . 
The interpretation of these experimental observations 
are provided by methods based on the RKKY interac- 
tion [7,8], Hubbard type model [9], s-f exchange model 
[10] and ab-initio total energy calculations [11-13]. Naive 
application of the RKKY method (assuming spherical 
Fermi surface and uniform distribution of spins in the 
ferromagnetic layers) leads to small oscillation periods 
[7]. Hence Bruno and Chappert [8] extended the gen- 
eral theory of the RKKY method and got satisfactory 
oscillation periods for Co/Cu/Co and Fe/Cu/Fe multi- 
layers considering the topological properties of the Fermi 
surface of the spacer layer and the moment distribution 
within the ferromagnetic layers. Using a Hubbard like 
one band model [9], Edwards et. al. showed that the 
exchange coupling exhibits long period oscillations with 
respect to the spacer layer thickness for certain posi- 
tions of the Fermi level. Urbaniak-Kucharczyk used the 
s-f exchange model to study the interlayer coupling ex- 
pressing the effective exchange integrals in terms of the 
electron susceptibility [10]. Convincing results for the 
fcc-(lll),(100) and (110) spacers are found. Ab-initio 
total energy methods are restricted to small spacer layer 
thickness as they become very difficult when the spacer 
thickness is large. An extensive discussion of the gen- 
eral theory of the interlayer exchange is given by Bruno 
[14] and the illustration of the theory for the case of 
Co/Cu/Co multilayer is also presented by him. The in- 
terlayer exchange coupling is found to have significant 
influence on the magnetic properties of the ferromagnetic 
sublayers. Ney et. al. found an oscillating behaviour in 
the Tc of Co/Cu/Ni trilayers with respect to the spacer 
layer thickness [1] and a theoretical analysis of this effect 



based on the Hubbard model is reported by Wu et.al. 
[15]. 

Our interest is to study the exchange coupling of the 
ferromagnetic (M) layers in M/N/M multilayers with M 
being a local moment ferromagnetic metal and N being 
a nonmagnetic metal. In the local moment metals, the 
magnetic properties are dominated by the intra-atomic 
exchange interaction acting between the conduction elec- 
trons and the local moments (s-f exchange) and the mag- 
netic coupling between the local moments situated at the 
lattice sites is mediated by the conduction electrons. The 
above mentioned correlation effects are properly taken 
into account in the s-f model which is often referred to 
as the ferromagnetic Kondo lattice model in the recent 
liturature and the modified RKKY method [16,17] de- 
scribes the exchange coupling of the localised spins me- 
diated by the correlated conduction electrons. Thus this 
approach leads to the self-consistent evaluation of the 
exchange integrals acting between the local moments sit- 
uated at the lattice sites and hence we use this s-f model 
along with the modified RKKY method to study the in- 
terlayer exchange coupling in M/N/M multilayers with 
M being a local moment metal. The starting Hamil- 
tonians of our method and the method of Urbaniak- 
Kucharczyk [10] are similar but the approaches are com- 
pletely different (see sect.V). Further, our method and 
the Urbaniak-Kucharczyk's method are different from 
the Bruno's method [14] in the sense that in the Bruno's 
method, the lEC is derived from the interference effects 
of the electron waves whereas in our method and the 
Urbaniak-Kucharczyk's method, the lEC is derived using 
the indirect coupling of the localised spins mediated by 
the conduction electrons. In our method the full Green 
function of the s — f system is evaluated and the inter- 
layer exchange coupling (lEC) is studied using the above 
said modified RKKY method. In the present work, we 
have carried out our study for a model system with only 
one conduction band per layer and we hope that this 
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study will form the basis for the calculation for a real 
system. The influence of the intcrlayer exchange cou- 
pling on the magnetic properties of the ferromagnetic 
sublayer can also be studied using this method. However 
this study is planned for a future paper. In this paper 
we present our results on the intcrlayer exchange cou- 
phng in the M/N/M multilayer at T=OK with different 
thicknesses of the non-magnetic spacer layer. We have 
further studied the influence of the conduction electron 
concentration on the lEC. We present our results in the 
following sections along with a brief discussion of the s-f 
model and the modified RKKY method. 



II. S-F MODEL FOR THE FERROMAGNETIC 
FILMS 

In this section we shall present a brief discussion of the 
s-f model. We shall present the theory for an 'n' layer 
film and the theory can then easily be adapted for the 
required M/N/M multilayer geometry by assuming that 
there are no localised spins in the non-magnetic spacer 
layers. The theory and the mathematical formulation of 
the s — / model are described in many of the earlier pub- 
lications [16,17]. Hence we will not present those details 
here. However, for the sake of completion, we will present 
a brief summary of the model. We will consider a ferro- 
magnetic film with n layers. The film is characterised 
by a two dimensional Bravais lattice vector having an 'n' 
atom basis. The 'n' atom basis corresponds to the 'n' 
layers of the film. A lattice vector of the film may be 
given as 

Ri is the two dimensional Bravais lattice vector and 
is the basis vector. At each site of the film a localised 
spin 'S' is present and the film is assumed to have only 
one conduction band per layer. The electrons in the con- 
duction band are exchange coupled to the local moments. 
This situation is well described by the s — f model, the 
model Hamiltonian of which is given as 



The first term describes the conduction electrons and the 
second term represents the interaction of the conduction 
electrons with the local moments. i,j represent the site 
indices of the two dimensional lattice and a, [3 represent 
the layer indices, ff^^ is the hopping integral and j is the 
s — f exchange coupling strength. All the information 
concerning the electronic structure and magnetic prop- 
erties of the system described by the above said Hamil- 
tonian can be obtained from the retarded single electron 
Green function 
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Evaluation of this Green function is proceeded by the 
equation of motion method. The equation of motion is 
written as 



Defining the selfenergy MijdE) as 

and making a Fourier transformation (with respect to the 
spatial variables) of the equation of motion leads to 
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In the matrix form the above equation may be written 

n 



as 



e(k) is the Bloch energy matrix and for an 'n' layer 
film it will take the form 



£(k) = 



V£"l(k) £"2(k) 



£l"(k)\ 

e2»(k) 
e""(k) ) 



The diagonal elements are the intralayer Bloch energies 
and the off-diagonal elements are the interlayer Bloch 
energies. 

Mcr{E) is the self-energy matrix and it takes the form 
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The diagonal elements are the intralayer self-energies and 
the off-diagonal elements are the interlayer self-energies. 

The layer dependent spectral density may now be ob- 
tained from the Green function matrix as 



S^^{E) = --ImGl^{E) 
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The layer dependent quasiparticle density of states 
(QDOS) is 



The evaluation of the self-energy is described in many 
of the earlier publications [18-20]. Once the self-energy 
is calculated, the matrix Green function of the film can 
immediately be obtained. 



III. MODIFIED RKKY METHOD 



As the exchange integrals are dependent on the Green 
function of the s — f system, it is obvious that they are 
dependent on the electron self-energy. 

It is possible to obtain the conventional RKKY interac- 
tion from this formalism through the first-order approx- 
imation 



{0)al3, 



In this case J"''(q) will be identical to the well-known 
RKKY expression. 



The modified RKKY method essentially aims at the 
evaluation of the indirect exchange coupling between two 
localised moments mediated by the conduction electrons 
which are coupled to the localised spins situated at the 
lattice sites through the s — / exchange acting between 
them and the localised spins. In order to obtain this ef- 
fective exchange interaction between the localised spins, 
the s — f interaction is mapped to an effective Heisen- 
berg Hamiltonian by averaging out the conduction elec- 
tron degrees of freedom [17]. This procedure allows us 
to introduce an effective exchange coupling between the 
localised moments as a functional of the conduction elec- 
tron self-energy. The details of the method may again be 
found in the earlier publications [17] and hence we will 
not elaborate the method here. However for the sake of 
completion, we will quote here the required results which 
are used in the calculations. The exchange integral acting 
between two localised moments is given by this method 
as 

q 

where J"'^(q) is given as 
and 

-g j^j^gg electron Green function and G'^^{E) is 
the already defined Green function of the s — / system. 
/_ {E) is the Fermi function. The q dependent exchange 
integrals, thus can also be written in a matrix form for 
an 'n' layer film as 

/J"(q) Jl2(q) .... jl„(q)x 

. J2i(q) J22(q) .... J2n(q) 

J(q) = 

Vj"i(q) J"2(q) .... J""(q)/ 



IV. INTERLAYER EXCHANGE COUPLING 

Our primary interest is to evaluate the indirect ex- 
change coupling between two local moment ferromagnetic 
layers separated by a nonmagnetic spacer layer including 
the effects of the above mentioned s — / exchange acting 
between the conduction electrons and the localised spins 
of the ferromagnetic layers. Hence we have made use of 
the modified RKKY method to evaluate the required ex- 
change integrals. The matrix Green function of the 'n' 
layer film can easily be converted to describe the M/N/M 
multilayer by assuming that there are no localised spins 
in the spacer layers and hence the Green functions of the 
spacer layers will be just the free electron Green func- 
tions. The matrix Green function of the M/N/M multi- 
layer thus can be generated and it can subsequently be 
used to obtain the exchange integral matrix J(q) dis- 
cussed in the previous section. The exchange interaction 
acting between two moments located in the bottom and 
topmost ferromagnetic sub-layers {Jij) of the multilayer 
can now be calculated from J(q). Summation over all the 
exchange integrals {Jij) acting between the localised mo- 
ments situated at the lattice sites of these two ferromag- 
netic layers will yield the interlayer exchange coupling. In 
all our calculations, we assumed the geometry of two fer- 
romagnetic monolayers separated by a spacer layer and 
the spacer layer thickness was varied upto thirty mono- 
layers. The geometry of the multilayer with a spacer 
thickness of five monolayers is shown in Fig.l. 

The Bravais lattice of each layer is assumed to be a 
squre lattice and the nearest neighbor intralayer and in- 
terlayer electron hopping is assumed to occur in the mul- 
tilayer. The magnitude of the hopping is assumed to be 
the same in all the layers. The value of the s — / coupling 
constant j is assumed to be 0.2 eV and the magnitude 
of the localised spin is taken as 7/2. The electron self- 
energy is actually dependent on the band occupation and 
temperature and hence the exchange integrals will also 
be dependent on the band occupation and temperature. 
We have studied the influence of the band occupation on 
the exchange integrals keeping the same band occupation 
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in all the layers and the influence of temperature on the 
interlayer exchange coupling is planned for a forthcoming 
paper. We have further evaluated the interlayer coupling 
using the conventional RKKY method and the results of 
our calculations are discussed in the following section. 




Fig.l Geometry of the M/N/M multilayer with a 
spacer thickness of 5 monolayers 



V. DISCUSSION OF RESULTS 

The calculated interlayer exchange coupling (lEC) us- 
ing the modified RKKY method at two different band 
occupations (n=0.2, 0.8) are presented in Figs. 2 and 
3. The lEC oscillates between ferro- and antiferromag- 
netic configurations. This behaviour is in agreement with 
that of the trend seen in many experimental works. As 
it is a model calculation, the oscillation periods and na- 
ture of the magnetic coupling with respect to the spacer 
layer thickness cannot be compared quantitatively with 
the experimental works. Hence our next goal is to ex- 
tend the calculations to a real material film where the 
calculational complexity however will be enormous. 
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Fig. 2 The interlayer exchange coupling (lEC) obtained 
from the modified RKKY method at various values of the 
spacer thickness at a band occupation of 0.2 (T=OK) 
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Fig. 3 The interlayer exchange coupling (lEC) obtained 
from the modified RKKY method at various values of the 
spacer thickness at a band occupation of 0.8 (T=OK) 

Our interest in this study is to point out the impor- 
tance of electron correlation effects (self-energy) in de- 
ciding the nature of the magnetic coupling of the ferro- 
magnetic layers with respect to the spacer layer thick- 
ness and in influencing the oscillation periods. Hence 
we carried out the calculations using the conventional 
RKKY method also. The results obtained using the con- 
ventional RKKY method for the same band occupations 
are presented in Figs. 4 and 5. It may be seen from 
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the figures that the influence of correlation on the mag- 
netic coupUng is significant. When the spacer consists of 
one monolayer, the coupling is antiferromagnetic for the 
uncorrelated film (Fig. 4) whereas it is ferromagnetic for 
the correlated film (Fig. 2). The behaviour of the mag- 
netic coupling at larger spacer layer thicknesses also show 
significant differences between the uncorrelated and cor- 
related films. The magnitude of the coupling strength 
is also drastically modified because of the correlation ef- 
fects and this may also be seen from the figures. Thus our 
studies clearly demonstrates the influence of the correla- 
tion effects on the exchange coupling of the ferromagnetic 
layers. 
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lEC as a function of the s — / exchange coupling strength 
(j) for the multilayer with a single spacer layer using 
the modified RKKY as well as the conventional RKKY 
methods. The results are plotted in Fig. 6. It may be 
seen from the figure that, both the methods predict an- 
tiferromagnetic exchange coupling at small values of j. 
However, after a critical value of j, the modified RKKY 
method predicts a transition of the exchange coupling to 
ferromagnetic nature, whereas the conventional RKKY 
method predicts that the antiferromagnetic coupling is 
retained for all values of j. The effect of correlation on 
the lEC is thus clearly demonstrated in Fig. 6. The mod- 
ified RKKY method further predicts that the lEC gets 
saturated after certain value of J, whereas in the case of 
conventional RKKY method, the lEC can never reach 
saturation. As the modified RKKY method uses the full 
Green function of the system which properly takes into 
account of the s — f interaction, it is obvious that these 
calculations are much more reliable than the conventional 
RKKY method. On the other hand, in the method of 
Urbaniak-Kucharczyk the effective exchange integrals are 
expressed in terms of the electron susceptibility and the 
s — f interaction strength j. The expression for the ef- 
fective exchange integrals is quadratic in j and hence the 
behaviour of the lEC seen in Fig. 6 is not abtainable from 
his method. The lEC is found to be sensitive to the band 
occupation also and it may be seen from the Figs. 2, 3 
and 4, 5. The band occupation also alters the nature and 
magnitude of the lEC. 



Fig. 4 The interlayer exchange coupling (lEC) obtained 
from the conventional RKKY method at various values of 
the spacer thickness at a band occupation of 0.2 (T=OK) 
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Fig. 6 Variation of the interlayer exchange coupling 
(lEC) with respect to the s — / exchange coupling 
strength (j) for the M/N/M multilayer with a monolayer 
spacer (T=OK) 



Fig. 5 The interlayer exchange coupling (lEC) obtained 
from the conventional RKKY method at various values of 
the spacer thickness at a band occupation of 0.8 (T=OK) 

In order to study the evolution of the correlation effects 
on the interlayer coupling, we have further calculated the 



In order to obtain the oscillation periods, we performed 
a Fourier analysis of our data. The dotted lines shown 
in Figs. 2-5 are the Fourier fit to our data. The Fourier 
analysis of our data corresponding to the modified RKKY 
method revealed that the oscillation period at n=0.2 is 
3.22 whereas as it is 4.83 when n=0.8. Thus we found 
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only a single period oscillation in the lEC. The oscillation 
periods are almost the same for the RKKY data also. It 
may be seen from Figs. 2-5 that there is a overall sign 
change of the decaying wave from the modified RKKY 
to conventional RKKY for n=0.2, whereas it is absent 
for n=0.8. The magnitude of the oscillation is very much 
reduced in the modified RKKY compared to that of the 
conventional RKKY. 



VI. CONCLUSIONS 

We have evaluated the interlayer exchange coupling of 
two local moment ferromagnetic sublayers in a M/N/M 
multilayer at different spacer thicknesses. The calcula- 
tions were carried out using the s — / model and the 
modified RKKY method which provide a self-consistent 
description of the exchange coupling between the mo- 



ments situated at the lattice sites of the two ferromag- 
netic sublayers. The lEC is found to have an oscillating 
behaviour with respect to the spacer layer thickness and 
it oscillates between ferro- and antiferromagnetic config- 
urations. Inorder to demonstrate the infiuence of correla- 
tion effects on the lEC, we have further evaluated it using 
the conventional RKKY method too. We find significant 
influences of the electron correlation effects on the lEC 
and the correlation effects are found to change the nature 
and magnitude of the exchange coupling. The variation 
of the lEC with respect to the band occupation is also 
studied using both the methods and the band occupation 
also found to alter the nature and magnitude of the lEC. 
The oscillation period of the lEC is also found to de- 
pend on the band occupation. Our method is very much 
sophisticated than the conventional RKKY method and 
hence a calculation for real local moment multilayers can 
be expected give reliable and interesting results. 
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